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1.0  INTRODUCTION 


This  is  the  fourth  and  final  quarterly  prog‘*  *,ss  report  for  a  one  year 
study  of  convergence  in  adaptive  radar.  During  *  series  of  earlier  adaptlv 

^  I 

radar  studies''  *  ■'  at  TSC,  funded  by  Naval  Air  stems  Command,  a  method 
of  optimizing  both  the  angular  and  doppler  res^.,  nsc  in  AMTI  (airborne 
moving  target  Indication  radars  was  described  •’Jialyzed,  and  simulated 
in  some  detail.  These  studies  showed  that  ai: ..  tive  AMTI  radar  can 
provide  a  major  Improvement  in  performance  '  c  a  mixed  clutter/ECM 
environment. 

In  the  first  progress  report  the  b,,sic  adaptive  AlITI  technique  was 

described,  and  three  adaptive  algorithms  were  compared.  These  algorithms 

were  sample  matrix  inversion  SMI,  inverse  matrix  update  IMU,  and  the 

Applebaum  loop  algorithiii.  Both  SMI  and  IMU  algorithms  are  shown  to  be 

very  fast  in  terras  of  number  of  data  samples  required  for  convergence. 

These  algoritlans  are  more  complicated  than  the  loop  technique,  though, 

and  are  probably  best  Implemented  using  digital  techniques.  Section.  2 

■of  Jthls  report  present^  results  of  a  computer  simulation  of  analogue  to 

digital  (A/D^  conversion.  A  theory  predicting  the  performance  loss  due 

"  . 

to  quantization  was  presented  in  progress  report  3  of  this  contract  and 
is  included  as  Appendix  I.  This  theory  is  compared  with  simulation  results 
in  Section  2.  The  program  used  to  obtain  the  results  of  Section  2  is  ^ 
included  as  Appendix  II.  In  Section  3  of  this  report  some  further  results 
are  shown  for  the  IMU  technique,  using  the  program  listed  in  reference  3. 
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2.0  SIMULATION  CF  QUANTIZATION  AND  LPIITIKG 

In  this  section  we  consider  the  effect  of  converting  the  Input 
analogue  data  to  digital  data  In  an  adaptive  AMTI  array,  Including  the 
ability  to  reject  clatter.  A/D  conversion  introduces  two  different 
types  of  non  linearity,  one  is  quantization  which  produces  a  step  like 
approximation  to  a  straight  line  transfer  function  and  the  other  is 
llmitirg  which  renders  the  output  constant  after  a  particular  input 
level  is  reached.  These  effects  are  important,  since  the  number  of 
bits  used  in  a  quantizer  is  limited  for  economic  reasons  and  must  be 
as  small  as  is  consistent  with  good  performance.  The  above  non- 
linearities  are  simulated  separately  and  together  using  the  computer 
programs  in  Appendix  II.  The  simulation  of  quantization  and  limiting 
separately  allows  comparison  with  the  theoretical  results  for  quantization 
noise  given  in  Appendix  I. 

The  example  considered  for  simulation  is  an  AMTI  radar  system  in 
which  the  weights  ere  multiplied  by  the  output  of  a  number  E  of 
receiving  antenna  elements  on  each  of  a  number  I  of  pulses  and  then 
added  to  obtain  the  radar  output.  The  non  linearities  are  introduced 
as  shown  in  Figure  2.1.  The  weights  as  in  earlier  studies  can  be 
nearly  optimized,  using  data  samples  collected  by  the  radar.  A  single 
sample  of  data  is  defined  as  the  set  of  returns  frco  a  single  range 
bin  at  each  element  on  each  pulse,  so  that  for  a  2  pulse  MTI  with  A 
antenna  elements,  c,  sample  would  contain  eight  complex  values.  All 
Wfights  and  voltages  are  complex  thereby  preserving  both  phase  and 


Figure  2.1.  Functional  Diagram  of  Adaptive  Radar  With  Non-Linearities. 
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amplitude  information.  Each  channel  shown  in  Figure  2.1  consists,  in 
fact,  of  two  quadrature  channels,  one  in-phase  1,  the  other  90®  out- 
of-phase  Q.  On  transmission  the  elements  are  combined  with  equal 
weighting  to  form  a  fixed  transmit  beam. 

In  a  computer  simulation  where  the  ii;terference  or  noise  (in  this 
case  distributed  ground  clutter  and  receiver  noise)  is  known  exactly, 
an  ideal  covariance  matrix  M  for  a  linear  system  is  computed  as 


J 

M  -  £ 
K,  jti 


- 


(2-1) 


where 

V  *  the  voltage  at  the  m*"^  element  (pulse)  from  the 
j  scatterer 

D  ■  distance  from  scatterer  to  element 

X  ■  transmitter  wavelength 
*  means  complex  conjugate 

For  convenience  the  indexing  over  antenna  elements  and  pulses  is  combined 

sc  that  V  becomes  V  where  m  “  e  +  (p-l)E.  Included  in  V  is  the 
ep  m  '.t'  / 

scatterer  amplitude  and  any  other  amplitude  weighting  such  as  the  trans¬ 
mitter  antenna  pattern  and  individual  element  patterns. 

Hie  signal-to-cl'itter  or  HTI  gain  it.  given  in  vector  notation  by 


G 


Mi?. . 

WMW 


(2-2) 
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for  a  set  of  weights  W,  where  (S/C)^  is  the  signal  to  clutter  obtained 
on  a  single  pulse  with  uniform  weighting  on  receive  and  transmit.  This 
factor  is  used  for  normalization. 

fl-41 

As  explained  in  previous  references  the  weights  that  maximize 

S/C  or  G  are  known  to  be 


-1  * 

W  -MS 
o 


(2-:) 


and 


IW  SI 
o 


"  WMW* 
o  o 


1  “1*1 
*=  SM  S 


(S/C), 


(2-4) 


where  S  is  the  desired  signal  vector.  When  non-linearities  are  present, 
for  example  as  indicated  in  Figure  2.1,  the  clutter  covariance  matrix 
cannot  always  be  readily  calculated.  Instead  it  can  be  found  by  using 
computer  generated  data  samples  to  form  a  sample  covariance  matrix  (or 
its  inverse)  as  it  would  be  done  in  an  adaptive  radar  using  digital 

A 

processing.  The  eample  covariance  matrix  M  is  given  by 


M 


mn 


1  A  * 

“  )  U  U 

S  jti  ni 


(2-5) 


where  the  summation  is  over  independent  samples  S.  The  Ujjj 
element  voltages  after  the  non-linearities  are  applied.  In  the  computer 
simulation,  each  voltage  for  a  given  sample  is  generated  by  choosing 
a  different  set  of  Gaussian  distributed  random  numbers  for  c  prescribed 
set  of  scatterers.  The  sample  covariance  matrix  is  an  approximation 

[5] 

tf'  the  ideal  matrix  for  a  set  of  scatterers.  As  it  was  shown  previojsly  , 


the  convergence  is  very  fast. 
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The  samj-le  covariance  matrix  can  be  used  to  compute  a  set  of  weights. 


Tnus 


(2-6) 


The  question  then  is  how  to  test  the  weights  W  generated  from  a  sample 
covariance  matrix.  It  is  meaningful  with  a  linear  ryjitem  to  test  the 

A 

Wg  using  the  ideal  covariance  matrix  since  M  -♦M  as  the  number  of  samples 
becomes  large.  With  non-linearities  present  this  I.j  not  so. 

Two  cases  are  of  practical  Interest.  In  case  1  a  sample  matrix  is 
found  and  used  to  create  optimum  weights  W^^  for  M^.  Then  these  weights 
are  tested  on  the  same  sample  data  used  to  obtain  the  weights.  Thus  Case  1 


Wi  » 


(2-7) 


WiMiW*  j 


(2-8) 


where  and  are  the  covariance  matrices  with  and  without  signal. 


Case  2  is  similar  to  the  above  except  that  the  sample  matrix  used  for 
testing  the  weights  is  made  up  of  different  samples.  Thus  in  Case  2 


(S/C), 


(2-9) 

(2- 10) 


In  this  formulation  the  signal  is  not  present  in  the  simulation  and  is 
assumed  to  be  unaffected  by  the  non-linearities.  Since  signal  loss  is 


quite  small  even  ..1th  hard  limiting  this 


a  reasonable  apjn  u.'- Imat  Ion . 
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Case  ?  ±s  easier  to  apply  in  practice.  For  Case  1  to  be 
implemented,  the  data  from  all  elements  is  stored  and  used  to  form  a 
covariance  matrix,  then  the  weights  derived  are  used  on  each  range 
cell  individually.  This  requires  more  storage  than  in  Case  2,  where 
the  weights  are  used  on  succeeding  range  bins  as  the  data  arrives. 

In  practice  the  weights  would  be  applied  to  each  range  bin  individ¬ 
ually  to  maximize  the  signal- to-clutter  gain.  For  simulation  purposes 
one  is  Interested  not  in  individual,  but  in  average  performance,  there¬ 
fore  one  can  apply  the  weights  to  the  already  summed  covariance  matrices 
or  M2  as  in  Equation  2.8  or  2.10.  This  is  equivalent  to  applying  the 
weight  to  each  sample  individually  and  averaging  the  n -suits. 

The  relations  discussed  so  far  contain  the  assumption  that  the 

signal  is  not  present  during  the  collection  of  data,  i.e.  that  no  signal 

is  present  in  the  covariance  matrix  used  to  determine  the  optimum  weights. 

This  situation  cannot  be  strictly  true  in  Case  1,  or  no  signals  would  be 

detected.  It  could  be  true  in  Case  2,  if  the  data  used  for  determining 

the  weights  were  obtained  over  a  region  known  to  contain  no  signals. 

[5] 

From  previous  s  ■Emulations  it  has  been  found  that  little  degradation 
in  performance  occurs  when  signals  are  present.  If  signals  are  small 
the  dominance  of  clutter  which  is  present  xn  many  range  bins  compared 
to  the  one  or  two  for  signals  determines  the  weights  as  if  the  signals 
were  not  present.  If  the  signal  xs  large  detection  is  not  a  problem 
anyway,  so  that  Case  1  as  we’l  as  Case  2  apply  ».o  realizu^le  cas-es. 

Two  types  of  non-linearities  are  explored  r-  tne  elsrul ;tion , 
quantization  which  o:^-urs  in  the  conversion  of  ^•inalog  data  to  digital 
as  indicated  in  Figure  2.2,  where  a  linear  inaction  is  replaced  by  a 
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Figure  2.2.  A/D  Converter  With  Equal  Step  Sizes 
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stair  step,  and  limiting  as  also  indicated  in  Figure  1.1,  The  limitiiig 
can  be  either  in  the  I  and  Q  channels  in  which  case  phase  distortion 
occurs,  or  amplitude  (envelope)  limiting  as  might  be  performed  at  IF 
before  the  quadrature  channels  are  determined.  In  the  latter  case 
the  amplitude  is  limited  without  causing  phase  distortion. 

2.1  QUANTIZATION 

In  Table  2.1  are  shown  results  for  quantization  noise  with  no 
limiting  for  various  values  of  step  size  A  divided  by  RMS  input  voltage  o. 
The  example  used  consists  of  a  4  element  antenna  with  .5X  spacing  between 
elements  and  2  pulse  MTI.  The  beam  is  broadside  to  the  array  and  is 
aimed  in  the  iirection  of  travel  of  the  radar  which  is  moving  at  a  rate 
of  .2X  per  pulse.  Thirty  scatterers  are  distributed  uniformly  over  180* 
to  simulate  ground  clutter.  Under  these  conditions  the  ideal  MTI  gain 
as  defined  earlier  is  53.8  dB. 

The  losses  shown  are  the  MTI  gain  values  obtained  with  quantization 

divided  by  the  ideal  performance  and  expressed  i.i  dB.  In  this  Table  24  or 

[11 

3N  ■■  samples  are  used  to  form  the  sample  covariance  matrix  which 

A 

is  used  to  generate  the  weights  These  weights  are  then  tested  on  M^^ 

for  Case  1  and  on  which  is  formed  from  an  equal  number  bur  different 
set  of  samples.  The  experiment  is  repeated  12  times  in  this  Table,  i.e. 

12  runs  and  the  results  of  the  ,.ans  are  averaged.  The  variance  as  shown 
i^  quite  small. 

Some  of  the  simulation  vuL.es  shown  foi  Case  1  are  negative  indicating 

^^^Where  N  is  the  number  of  degrees  of  freedom  which  ?n  this  case  is 
2  X  4  8. 
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Table  2.1 

QUANTIZATION  LOSS  (dBj  (NO  LIMITING) 


2  pulse  MTI  4  elements 


30  clutter  scatterers.  uniformly 
distributed  over  180" 


Platform  motion  0.2X  per  pulse 

Beam  angle  0° 

from  velocity  vector 

Ideal  performance  (no  quantization) 

G  •=  53.8  dB 

Samples  used  in  M  *  24  (3N) 

Runs  averaged 

»  12 

I 

1 

Theoretical  Quant.  Loss  »■ 

Theoretical  Sample  Loss 

L  «=  1.43  dB 
s 

A/o 

1 

Case  Lq  Lq 

V^s 

Simulation 

Loss  Variance 

I 

0 

1 

-1.65 

.034 

2  0  0 

1.4 

1.24 

.021 

.01 

1 

.02 

.073 

r 

2  1.90  3.47 

3.33 

2.97 

.042 

.02 

1 

3.44 

.056 

2  5.05 

6.43 

•  6.586 

.038 

i’ 

.04 

1 

7.873 

.075 

2  9.91 

11.33 

10.94 

.126 

.08 

1 

13.65 

.059 

» 

2  15.54 

16.98 

16.63 

.060 

Case 

1-W^  tested  on  Case  2 

tested 

on  All  values  are 

in  dB  (-)  indicates 

» 

a  gain.  A/o  »  step  size/RMS  voltage, 


t 
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an  Improvement  over  the  so  called  ideal  linear  case.  This  is 
possible  when  the  weights  generated  from  a  particular  set  of  samples 
are  used  on  the  same  set  as  in  Case  1.  The  same  weights  give  poorer 
results  when  used  on  other  sets  of  data  generated  by  the  same  process. 

Two  theoretical  results  have  been  developed  for  loss  from 
ideal  performance.  One  loss  is  due  to  finite  sample  size  and  is  given 


-  10-  l«8io 


/  S-f2-N\ 

\  S+1  ) 


(2-11) 


where  S  is  the  number  of  samples  used,  and  N  is  the  degrees  of  freedom, 
i.e.  E  •  P  in  this  case. 

The  other  loss  due  to  quantization  was  derived  in  Progress  Report  #3 
and  is  reproduced  in  Appendix  1.  This  is  given  by 


*  -1 

1  T  SM  S 

Lq  -  10.  Log^Q  *  , 

S  M  Q. 


(2-12) 


where 


(2-13) 


Using  another  approximation,  a  second  formula  for  the  loss  was 
obtained  as 


Lq  ■=  10.  Log^Q  |l  -  ~ 


6  *  -1 
SMS 


(2-14) 


In  Table  2.1  the  simulation  and  theoretical  losses  are  shovm. 

It  Is  seen  that  there  Is  good  agreement  between  Case  2  simulation 
results  and  the  total  losses  based  on  sample  size  and  quantization 

t 

loss.  The  approximate  results  for  loss  do  not  give  such  good 
agreement . 

Table  2.2  shows  the  losses  for  various  numbers  of  samples  in  the 
covariance  matrices.  Theoretical  and  simulation  results  are  in  very 
good  agreement  and  the  variance  seems  reasonable  for  the  values  chosen. 

2.2  LIMITING 

In  this  section  limiting  without  quartization  is  explored  by  simula¬ 
tion.  Two  kinds  of  limiting  are  considered.  Amplitude  or  envelope  in 
which  the  vector  amplitude  is  limited  but  phase  is  not  distorted,  and 
limiting  in  the  IQ  channels  where  both  phase  and  amplitude  are  affected. 
These  two  types  are  illustrated  in  Figure  2.3. 

While  there  are  theoretical  results  for  hard  limiting,  it  is  difficult 
to  obtain  theoretical  results  for  linear  operation  with  the  various 
degrees  of  soft  limiting  used  in  this  case.  Simulation  results  are  shown 
in  Table  2.3  for  limiting  in  the  IQ  channels  and  in  envelope,  as  a  function 
of  limit  levels  divided  by  RMS  input  voltage  a. 

It  is  Interesting  to  note  that  for  strong  limiting  Csmall  values 
of  S/o),  IQ  limiting  seems  to  produce  the  greater  loss.  This  perhaps 
is  due  to  the  greater  phase  distortion.  At  values  of  S/o  producing  less 
limiting,  losses  seem  to  be  higher  when  aiaplitude  limiting  is  used. 

Also,  the  variance  is  quite  high  for  the  lower  limiting  values,  doth 
of  these  results  may  be  due  to  the  fact  that  for  larger  values  of  S/o 
theie  lb  less  limiting  in  IQ,  thereby  causing  less  chance  of  a  degraded 
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Table  2.2 

QUANTIZATION  LOSS  (dB)  (NO  LIMITING) 


2  pulse  MTI  4  elements 

Platform  motion  0.2X  per  pulse 
Ideal  performance  (no  quantization) 
Samples  used  in  M  *  24  (3N) 

I 

Theoretical  Quant.  Loss  •= 


30  clutter  scatterers  uniformly 
distributed  over  180“ 

Beam  angle  0“  from  velocity  vector 

S/C  -  53.8  dB 

Runs  averaged  “  12 

Theoretical  Sample  Loss  L 

6 


Sample=64(8N) 

Sample=32(4N) 

Sample“16(2N) 

Runs«4 

Runs =9 

Run3=19 

Lg=.5  dB 

Lg'=1.03  dB 

Lg=2.3  dB 

1 

Sim. 

A/o  Case  L^  L^ 

L-+L„  Loss  Var. 

_ _ 

1 

-.61  .007 

-1.275  .015 

-3.02  .068 

.001 

2  .024  .024 

.52  .33  .001 

1.06  .923  .026 

2.33  2.35  .073 

1 

1.22  .03 

.55  .025 

-1.1  .121 

.01 

2  1.91 

2.4  2.3  .015 

2.94  2.7  .035 

4.2  4.01  .094 

1 

16.9  .007 

16.13  .047 

14.35  .273 

-1 

2  17.42 

.  . . 

17.9  18.02  .012 

18.46  18.4  .028 

19.72  19.37  .102 

15 


result.  When  limiting  does  occur  considerable  degradation  results. 

This  is  usually  Indicated  by  drastically  increased  variance. 

The  percent  limiting,  shown  in  Table  1^3,  is  obtained  by  counting 
the  number  of  times  that  limiting  occurs  during  simulation.  The  reason 
that  tb''-'e  is  more  limiting  in  the  er.velope  case  (the  circle  in  Figure 
2.3)  is  probably  that  the  sample  vectors  have  more  '‘non-limit”  area  in 
the  complex  plan  when  IQ  limiting  is  used  (the  squaie  in  Figure  2.3). 

2.3  LIMITING  AKD  QUANTIZATION 

Slmulatior  MTI  gains  are  shown  in  Figure  2.4  a.id  2,5  for  quantization 
and  limiting  combined  as  a  function  of  normalized  saturation  level  S/o 
(a  =  RMS  input  voltage).  The  results  are  shown  for  8,  9  and  10  bits. 

As  S/o  increases,  limiting  decreases  but  the  quantization  step  size  A 
increases  so  that  there  are  optimum  values  of  S/o.  Since  quantization 
noise  degrades  performance  slowly  compared  to  limiting,  there  is  a  fairly 
broad  region  of  S/o  values  over  which  to  operate. 

There  is  little  difference  as  indicated  by  Figure  2.4  and  2.5 
between  IQ  and  envelope  limiting  in  the  region  of  small  limiting  as 
might  be  expected.  With  strong  limiting  envelope  limiting  shows  better 
but  highly  degraded  performance. 

In  Figures  2.6  and  2.7  results  are  shown  for  a  case  with  higher  MTI 
gain.  In  this  case  the  scan  angle  is  90®  to  the  flight  path,  and  the 
interpulse  motion  is  close  to  X/4,  the  velocity  and  angle  for  which 
perfect  platform  motion  compensation  is  possible.  MTI  gain  is  plotted 
in  this  section  rather  than  losses  to  illustrate  performance  at  dltiereut 
gain  levels.  The  Case  1  performance  iu  wliich  is  used  with  rather 
than  (Case  2),  exceeds  the  steady  state  or  ideal  performaace.  This 


Table  2.3 


Simulation  Results  for  Limiting  Only 
2  pulse  MTl  4  elements  30  scatters 

Platform  motion  *=  0.2X /pulse  Beam  angle  0.0 

Ideal  G  «  53.8  dB 
Sample  “  64  (8N)  Runs  ®  4 


I  and 

1  Envelope 

S/a 

Case 

Loss(dB) 

Variance 

X 

Limiting 

Loss(db)  Variance 

% 

Limiting 

1.5 

1 

30.5 

.015 

27 

21.3 

.039 

32 

2 

31.7 

.001 

11. k 

.014 

2.0 

1 

24.0 

.07 

9 

18.4 

.056 

14 

2 

26.7 

.17 

20.6 

.015 

2.5 

1 

12.3 

•'  .533 

2.66 

12.2 

.162 

4.33 

2 

15.9 

1.25 

16.1 
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is  possible  since  the  weights  are  detenained  for  a  particular  set  of 
samples  and  used  on  the  same  sample  set. 

In  Figure  2.8  the  same  example  is  used  as  in  Figures  2.6  and  2.7 
but  the  range  of  S/a  is  extended  and  different  numbers  of  bits  are 
used.  This  Figure  shows  clearly,  that  while  there  is  an  optimum 
ratio  of  limit  level  S  to  average  voltage  a,  the  optimum  is  quite 
broad,  and  performance  declines  slowly  with  increasing  values  of  S/a. 
This  shows  that  it  is  clearly  better  to  operate  in  a  region  of  low 
limiting,  and  to  bias  the  operating  point  toward  the.  high  side  of  S/a 
when  a  is  changing  rapidly  or  is  not  accurately  known. 

Figure  2.9  shows  similar  results  for  a  lower  performance  level. 


4  elementB  2  pulses 

Element  spacing  ■  .5X  Intsrpulse  notion  • 

Samples  32  (AN)  Suns  «  4 

Scan  Angle  ■  90* 
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3.0  SOME  CONVERGENCE  EXPERIMENTS 
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In  the  first  progress  report  the  convergence  rates  of  three 

algorithms  were  compared  using  a  computer  simulation.  The  sample  matrix 

inversion  (SMI)  and  inverse  matrix  update  (I'MU)  techniques  were  shown 

to  be  comparable  and  very  fast  when  measured  by  t''c  number  of  independent 

data  samples  required  to  converge.  The  Applebaum  loop,  .jhile  much  simpler 

to  implement,  converges  much  more  slowly  requiring  hundreds  of  samples,  in 

the  example  shown,  compared  to  two  to  four  times  the  number  of  degrees  of 

freedom,  N,  required  by  SMI  and  IMU. 

The  SMI  technique  is  known  to  require  the  least  number  simples 

under  all  conditions,  and  requires  only  N(N+l)/2  complex  multiplications 
^'[5] 

to  form  M  \  These  properties  are  invariant  over  all  problems,  thereby 
making  the  technique  extremely  powerful. 

A  2 

To  invert  tht  matrix  M  for  SMI  requires  about  N  /3  multiplications 
or  divisions.  There  are  several  ways  of  doing  this,  i.e.  solving  the  set 
of  siiiiultaneous  equations  required  to  obtain  the  weights,  or  inverting 
the  matrix  by  ■''arious  methods.  Computational  problems  may  be  involved 
in  the  SMI  or  IMU  technique  when  high  speed  and  limited  numbers  of  bits 
is  a  requirement. 

The  DIU  algorithm  is  simple  to  program,  requiring  only  about 
2 

1.75  N  +  2.75  N  complex  multiplications,  per  data  sample.  This  can  be 
seen  from  the  recursive  equation 
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rl  1  ,-l  -c 

j+1  (1-a)  j  1-a  (i_a)+a(V^l“lv*) 


(3-1) 


where  the  cample  matrix  Is  given  by 


+  “''/jT 


(3-2) 


-1  *  2  , 
Basically  V  Is  formed  using  N  multiplications,  Is  formed 

by  taking  Its  conjugate,  and  to  obtain  ^  requires  a  further  N  multi- 

“1  —1  *  —1 
plications.  Since  is  Hermitian,  multiplying  V  by  requires 

N(N+l)/2  multiplications.  Finally  multiplying  by  the  constants  provides 
for  the  reamining  multiplications. 

The  simplicity  of  implementing  IMU  makes  it  interesting  to  see  how 
many  samples  are  required  for  convergence  both  as  a  function  of  a  and  of 
the  starting  values  of  For  fixed  values  of  a  an  exponential  weight¬ 

ing  is  applied  to  the  incoming  data  samples  as  they  are  used  to  update 

the  matrix.  Wlien  a  =  Try*  equal  weighting  is  applied  to  the  data  samples 

J+I 

used  to  update  M 

The  starting  values  of  M  ^  are  less  easily  parameterized  since  they 
could  be  anything  from  ai  arbitrary  choice  to  a  choice  which  is  close  to 
the  correct  value.  To  explore  IMU  as  a  function  of  the  initial  matrix 
one  wishes  to  construct  a  problem  related  to  a  practical  radar  situa¬ 
tion  in  which  external  interference  (clutter,  jamming  etc.)  is  changing. 
Also  one  would  like  the  weights  to  change  In  such  a  manner  that  a  required 
perl oi-manct*  is  maintained.  For  example,  one  might  require  the  weights  to 


2o 


adjust  as  the  antenna  rotates  between  pulse  groups.  This  was  treated 

rsi 

in  the  first  progress  report  .  For  another  example,  one  might  require 

the  weights  to  change  with  the  changes  in  clutter  with  range  as  treated 
in  a  previous  TSC  study The  latter  might  be  caused  by  terrain  changes 
such  as  a  land  to  sea  Interface. 


In  these  experiments  antenna  step  scanning  was  assumed  in  such  a 

manner  that  the  weights  are  optimized  at  one  angle,  then  the  antenna  is 

rotated  to  a  new -angle  and  held  fixed  while  new  samples  are  taken.  In 

Figure  3.1  the  TMU  process  is  started  with  the  diagonal  having  the  same 

values  as  the  diagonal  of  the  inverted  matrix  for  the  previous  angle. 

The  inverse  matrix  is  then  updated  with  the  new  data,  using  a  »=  7^  , 

j+i 

i.e.,  giving  equrl  weight  to  each  new  sample. 

A  large  rotation  angle  5®  was  chosen  to  accentuate  the  missmatch. 

The  results  for  the  IMU  and  SMI  techniques  are  seen  to  be  quite  similar. 
The  loop  convergence  performance  is  also  shown  for  the  case  where  steady 
state  performance  was  reached  at  the  previous  angle.  The  gain/TAU  ratio 
is  adjusted  so  that  control  loop  noise  in  the  steady  state  increases  the 
total  output  noise  by  25%^^^. 

It  seemed  reasonable  that  if  the  entire  inverse  matrix  for  a  previous 
starting  angle  were  used,  that  the  results  would  be  equally  good  if  not 
better.  Figure  3.2  shows  that  this  is  not  the  case,  at  least  for  this 


example.  The  IMU  performance  was  considerably  degraded  when  the  ideal 
M^^  for  the  previous  angle  was  used  as  the  initial  matrix. 

Figures  3.3  and  3.4  show  that  using  the  previous  inverse  matrix 


diagonal  even  when  it  is  multiplied  by  .1  or  10  gives  extremely  good 
agreement  between  the  SHI  and  IMU  techniques. 
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In  Figure  3.5  -  3.7  different  values  of  a  are  shown  using  the 
diagonal  of  the  previous  updated  matrix. 

These  figures  indicate  that  constant  values  of  a  are  much  less 
desirable  than  the  varying  value  which  gives  equal  weighting  to  the 
samples.  This  is  what  would  be  expected  in  the  step  scan  example, 
assumed  in  this  case.  If  a  steady  scan  rate  were  used  and  data  were 
taken  from  many  range  traces,  then  an  optimum  constant  value  for  a 
cou?.d  be  found. 

The  simulation  program  ADAPTM6  listed  in  the  first  progress  report 
was  used  in  these  experiments. 
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APPENDIX  I 


LOSSES  IN  ADAPTIVE  RADAR  PEPvFORJlANCE  DJE  TO  QUANTIZATION 


PRELIMINARIES 

Analyses  of  adaptive  radar  performance  usually  neglect  the  detailed 
effp-'ts  of  quantization  noise.  If  estimates  of  degradation  in  perfor¬ 
mance  are  made,  they  often  depend  on  costly  computer  simulations., 

This  situation  is  improved  here  for  most  cases  of  practical  interest 
by  developing  formulas  for  computing  the  signal-to-noise  loss  in  adpative 
radar  performance  due  to  quantization. 

Assume  the  adaptive  array-radar  has  a  receiving  array  of  M  identical 
elements.  Let  2j^(t)  be  the  complex-valued  process  received  by  the  k-th 
element  for  k=l,2,...,M  .  In  a  pulse-sampling  radar,  pulses  are  transmitted 
periodically  with  pulse  repetition  period  T.  For  this  case  the  '  fleeted 

CT 

signals  from  an  object  at  fixed  range  P  =  —  are  proportional  to  ..he 

coefficient  of  reflectivity  of  that  object.  Here,  of  course,  T  is  the 

round-trip  time  required  for  a  pulse  to  travel  from  the  radar  antenna 

to  the  object  and  return,  and  c  is  the  velocity  of  light.  Such  a  radar 

samples  the  reflectivity  of  the  object  sequentially  at  times  t^,t2, . . . 

where  t^  =  .  T.  The  sampled  data  set  associated  with  these  sampling 

times  is  the  sampled  data  set  {z,  (t  )}  of  signals  received  from  an  object 

K  n 

at  range  R.  It  is  convenient  for  mathematical  purposes  to  represent 
this  sampled-data  set  as  the  column  vector  (matrix). 
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where  N  =  MxL  is  the  number  of  dimensions  of  the  vector. 

Radar  detection  involves  a  choice  between  two  hypotheses,  the 
noise-only  hypothesis  Hq  and  the  signal-plus-noise  hypothesis 
Assuming  additive  Gaussian  noise  of  zero  mean,  the  expected  value  of 
X,  given  hypothesis  Hq  (noise  only),  is 

EX  =  0. 

Similarly  the  expected  value  of  X,  given  (signal  plus  noise)  is 

EX  <=  S 


where  S  is  the  column  vector  (matrix) 
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of  signal  echoes  that  night  be  expected  from  an  object  at  range  R  at 
the  M  elements  of  the  array  from  the  L  pulses.  For  hypothesis 
the  noise  vector  is  given  by  N  =  X  -  S.  Assume  that  all  components  of 
M  are  Gaussian  and  jointly  distributed. 

In  order  to  detect  signal  vector  S  in  the  received  sampled-data 
vector  X  one  nust  design  a  filter  (a  linear  functional)  f-^r  vector  X 
which  is  tuned  to  signal  S.  Such  a  sampled-data  filter  is  the  scalar 

N 

y "  E  ”k  ’‘k ' 

k=l 


where  W  is  the  weight  vector 


W  <= 


of  cotaplex  numbers,  denotes  complex  conjugation  and  W  denotes  the 
conjugate  transpose  of  the  column  matrix  W. 


2.2  MOMENTS  OF  QUANTIZED  VIDEO 

To  study  the  effects  of  quantization  on  adaptive  radar  performance 

one  needs  to  find  the  moments  of  the  sample  data  filter  y,  given  in  (i). 

To  treat  this  problem  let  [x  ]  -  [u  )  +  ifv  ]  be  the  quantized  or  digi- 

n  n  u 

talized  value  of  the  complex  number  x  ==  u  +  iv  wheie  [v  ]  denotes 

n  n  n  n 

tlie  digital  value  of  the  am.  ginary  part  of  x^^.  In  the  conversion  from 


analogue  to  digital  let 
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[x  ]  •=  X  +  c(x  )  (2) 

n  n  n 

where 

e(x  )  *=  e(u  )  +  1  e(v  )  (3) 

n  n  n 


is  the  complex  error.  Here  e(u  )  is  the  real  part  of  e(x  )  and  e(v  ) 

n  n  n 

is  the  imaginary  part  of  error  e(x  )  in  the  A-D  conversion  process. 

n 

The  error  functions  eCu)  and  e(v)  are  sawtooth  functions  of  the  real 
and  imaginary  parts  of  x  =  u  +  iv,  respectively.  Explicitly,  e(u) 
and  e(v)  have  the  same  Fourier  series  representation,  given  by 


00 


k=l 


respectively,  where  A  is  the  quantizer  step  size  in  voltage  [see  Reference  1], 
For  signal-plus-noise  the  first  moment  of  [x^] ,  the  quantized 
version  of  the  n-th  component  of  the  sampled-data  vector  X,  is  by  (2) 
and  (3) 


E(x  ]  =•  e(x  +  e.(x  )) 
n  \  n  n  / 

“Ex  +  E(cx  ) 
n  n 

“8  +  Ee(a  )  +  1  Ee(v  ) 

n  n  u 


40 


Replacing  e(o  )  and  e(T^)  by  their  Fouriei  series,  eq.  (4),  and 

inverting  the  order  of  the  summation  and  expected  value  operator  E, 
yields 


E[xJ 


s  +  A 

n 


E 

k=l 


4  4  27rk  , 

i  sin  -7—  V  } 
^  n 


Since  and  are  independent  identically  distributed  Gaussian 

variatiea  with  means  RCs  )  and  I(s  ),  respectively,  where  R(s  )  and 

n  n  n 

l(s  )  are  the  real  and  imaginary  components  of  signal  component  s  , 
n  n 


E  sin 


1 


00 


“(u  -Rs  ) 
n  n 


e(u  )du 
n  n 


-2iiV(oM)^  /2lll«(s)\ 


and 


Taus 


X 


2TlkR(s  ) 

_ 

A 


+  i  sin 


2T!kI(s  )  \ 

__n_j 
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where  R(s^)  and  are  the  real  and  imaginary  components  of  the  n-th 

component  of  the  signal  vector  S  and  o  is  the  standard  deviation  of 
the  noise. 

For  most  cases  of  practical  interest  in  adaptive  radar  A  «  a. 


For  this  case 


E[x  1  *=  s  +  0 
n  n 


where  0(y)  denotes  the  order  of  y  as  y  tends  to  zero.  Evidently 
Elx^]  is  closely  approximated  on  the  average  by  the  s^,  i.e.  E[x^)  ~  s^. 

The  complex  covariance  matrix  of  [X],  the  quantized  version  of  the 
data  vector  X  for  noise  alone  is  by  (3) 

Mq«E[X][X]* 

«=  E(X  +  e(X))(X  +  E(X))* 

•=  EXX*  +  Ee(X)  {e(X)}*  +  E({e(X)}X*  +  X{£(X)}*) 

(5) 


where  *  denotes  conjugate  transpose  and  e(X)  is  the  column  vector  of 
complex  quantization  errors  e(x^),  as  given  by  (3),  namely. 


e(X)  ^ 
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The  first  term  of  (5)  is  the  covariance  matrix  M  of  the  sample-data 
vector  X,  i. e. , 


M  •=  rx 


To  treat  the  second  term  of  (5)  and  make  use  of  previous  results, 
we  must  transform  the  data  sample  vector. 


“2  +  ^''2 


“n  ^ 


«=  V  +  IV 


of  n  complex  components  into  the  data  sample  vector 


ot  2N  real  components  where  Is  the  real  part  of  sample  Is  the 

imaginary  part  of  sample  x^,  and 


and  V 


are  the  real  and  imaginary  parts  of  X,  respectively. 

Let  also  the  complex  covariance  matrix  M  and  its  inverse  M  ^  be 
decomposed  into  real  and  Imaginary  parts  as 


M^A+iB  and  M  =C  +  iD 


Then  the  2N  x  2N  symmetric  matrix 


A,  -  B 

B,  A 


has  the  inverse 


/C,  -  D 


where  Y  denotes  the  transpose  of  vector  Y. 


*  -1 

Next,  the  quadratic  form  X  M  X  can  be  shown  to  be  Identical  to  the 
T  -L  -1 

quadratic  form  Y  X  y  where  X  is  given  by  (9).  To  show  this 

X*m“Sc  *=  (U  +  iV)  (C  +  ID)  (U  +  iV) 

-  (U  -  iV)'^(C  +  iD)  (U  +  iV) 

=  (U  -  iv"^)  (C  +  iD)  (U  +  iV) 

T  T  T  T 

«=  U  CU  +  V  DU  -  U  DV  +  V  CV 

-1  T 

where  the  imaginary  part  vanishes  since  M  “  is  Hermitian  and  where  U 
denotes  the  transpose,  of  the  real  vector  U.  Similarly 

«=  (u’^C  +  v'^D,  -  u'^D  +  V^C)  j 

■=  u'^CU  +  V^DU  -  u’^DV  +  v’^CV 

Since  this  agrees  with  the  above,  the  identity 

X*m’'^  ■=  (11) 

is  true.  The  joint  probability  density  of  the  components  of  the  vector 
X  is 

P(X)  -  (Tt)“"lMr^  exp  (-X*m"^X) 


(12) 
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where  [m]  is  the  determinant  of  the  covariance  matrix  M.  Using  (11), 
this  density  is  easily  shown  to  be  equivalent  to  the  density 

P(Y)  »  (•.T)"^|Kl"^exp  (-yV4)  (13) 

where  Y  is  the  2n  component  vector,  given  by  (7)  and  K  ^  is  defined 
by  (9). 

The  above  remarks  will  now  be  used  to  evaluate  the  second  term  of 
(5).  The  (m,n)  element  of  this  is 


«=  E  e(u^)e(u  )  +  E  e(v  )e(v  ) 
m  n  m  n 

+  i(E  e(u  )e(v  )  -  E  e(u  )e(v  )^ 

\  m  n  n  n  / 

2(e  e(u  )e(u  )  +  1  E  e(u  )e(v  ))  (14) 

\  m  n  m  n  / 


where  the  last  line  follows  from  (9).  The  first  term  of  (14)  is 


E  e(u  )e(u  )  = 
m  n 


03  00 

Ji  "  “ 


)e(u  )e(u  )du  du 
m  n  m  n 


00  00 


k=l  £«1 


where 


Ikl 


// 


-  (u  +  u"  -  2p  u  u  )/2  (1-p  ) 

m  n  cm  m  n  /  mn  . 

e  >  sin 


/2ck  > 

U 


+  sin 


(^“n) 


du  du 
m  n 


J 1-P^ 


2iio  /  1-p 


cm 


a  E(u  u  )=  E(v  V  ) 

cm  m  n  m  n 


and 


cm 


P  =  1. 
nn 


The  above  integral  can  be  evaluated  to  yield  finally 


E{e(u  )e(u  )} 
la  n 


Z 


iill 


K+)l 


U 


K 


sin{4ir^ 


mn 


Again  for  the  cases  of  interest  A  «  a  so  that  to  first  order  in 


2  2 

e  only  the  k>=)l  terms  need  be  kept.  Thus 


00  2 

vi  f  \  f  \  V  1  (t)  ) 

E{e(u^)e(u^)  ^  ^  ^  \^J 

K=1  ^ 


(16) 


^  ,  -2TT^a/A)2 
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Finally  if  —  »  quantization  noise  power,  is  set  below  the  natural 
receiver  noise  of  the  radar,  only  the  terms  in  (16)  for  which  m=n  will 
persist.  To  this  final  approximation 


E4^o(, 


-2-n^  (a/A)^^ 


Ee(u^)e(u^)  =  0  + 


-2Tr^(CT/A)^(l- 


*^iDn^^  , 


using  the.  well-known  identity. 


1  u'- 


If  the  quantization  noise  power  is  again  set  helow  the  receiver 
noise,  the  imaginary  part  can  be  estimated  in  a  similar  manner  to  be 


Ee(u  )e(v  )  -  0  + 
n  n 


for  all  m  and  n  where 


a  y  -Euv  =-Evu 
mn  m  n  n  m 
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In  this  case  P  -  0,  since  by  (7),  (8)  and  (9), 
mn 


T  T  T 

B  “  EV  =  -EUV  =  “B 

so  that  B  Is  a  skew  symmetric  matrix.  Combining  ^17)  and  (18)  yields 

*  tP" 

E  E(X){e(X)}  *= -g-  I  (19) 


where  I  is  the  N  x  N  identity  matrix  as  an  estimate  of  the  seconr,,  term 
of  (5),  assuming  the  quantization  noise  is  comparable  to  the  receiver 
noise  power. 

A  similar  analysis  will  show  (see  reference  1)  that  tne  matrix 
elements  of  the  third  term  of  (5)  are  zero  to  the  order  cf  e 
Combining  this  with  (19),  yields  finally, 


-2iT^(a/A)^ 


a2 


Mq«=E[X][X) 


(20) 


2  2 

to  the  order  of  e**^*^  (1  This  result  will  now  be  used  to 

estimate  losses  in  detection  sensitivity  an  adaptive  radar  will 
suffer  as  a  function  of  the  quantization  step  A.  Better  approximations 
to  Mq,  then  given  by  (20),  particularly  for  a  larger  quantization  step  A, 
will  be  the  subject  of  a  future  study. 


SIGNAL- TO-NOISE  LOSS  DUE  TO  QUANTIZATION 

A  best  sampled-data  filter  of  form  (1)  is  one  in  which  the  weight 
vector  W  is  chosen  to  maximize  the  signal-to-r.oise  (S/N)  ratio.  The 
weight  vector  which  achieves  the  optimum  S/N  ratio  is  well  known  (see 
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for  example  Reference  2)  and  given  by 

W  =  m"^S  (21) 

o 

where  M  is  the  covariance  matrix  for  noise  only.  Since  M  is  often  not 

A 

known  a  priori,  an  estimate  M  is  used  Instead,  depending  on  K  samples  of  the 
noise  process.  One  estimate  is  the  sample  average.  If 
are  K  independent  samples  of  the  noise  process,  then  the  sample  average 
estimate  of  M  is  given  by 

^  "  K  £  (22) 

j=l 

If  M  is  used  in  (21),  then 

W  «  M"^S  C23) 

is  a  near  optimum  set  of  weights  to  be  used  in  filter  (1),  The  output 
of  this  filter  is 

A.  /V* 

u  *  W  X  (24) 

A 

where  W  is  given  by  (23).  The  output  signal-to-noise  ratio  conditioned 

A 

on  a  knowledge  of  W,  is 

(s/nIw)^  = 


Var(5T) 


*■^-1  . 2  .  1  '^-l 

(S  M  s:  /S  M  -^1  S 


(25) 
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I 


Previously  it  was  shown  [2]  that  if  this  signal-to-noise  ratio  was 
normalized  with  respect  to  its  maximum  value  that  this  quantity, 

namely, 

P(M)  =  (S/n1w)^  sV^S  (26) 

was  a  random  variable  which  was  in  the  interval  0  <  p(M)  <  1.  It  was 
further  found  that  this  normalized  signal-to-noise  ratio  has  a  probability 
density  P(P)  which  depedned  only  on  N,  the  number  of  components  of  X, 
and  the  munber  K  of  sample  vectors. 

A 

Previously  in  Reference  2  the  quantity  p(M)  was  the  (S/H)  loss  ratio 

A 

which  depended  on  M,  the  estimate  of  M,  without  any  assumption  of 

A 

quantization.  To  make  this  ratio  depend  on  quantization  let  be  the 
sample  average  estimate  of  M,  including  quantization,  i.e.  let 

E  (27) 

3=1 

where  denotes  the  quantized  value  of  sample  vector  Then 

the  normalized  signal-to-noise  ratio  for  this  estimate  of  M  is 


^Total 


(Mq) 


ss 


sV^s 


("'"I  Vo  \ 

s‘Hp^S  \sV^S  / 
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-  pcSq) 


] 

\sV^s/ 


(28) 


where 

Mq  =  E[X][X]* 

is  the  complex  covariance  matrix  of  [X],  given  by  (5)  and  (20).  For 
K  reasonably  large  the  statistics  of  [X)  will  be  close  to  Gaussian, 
hence  the  statistics  of  P(Mq)  will  closely  approximate  the  statistics 
of  p(M)  as  given  by  (26).  Thus,  taking  the  expected  value  of  (28), 
yields  by  equation  (18)  of  Reference  (2), 


^^Total 


\  •=+!  \s%-h  j 


(29) 


This  is  the  expected  loss  in  S/N  ratio,  due  firstly  to  the  fact  that 
only  a  finite  number  K  of  samples  were  used  to  estimate  M,  and  secondly 
to  the  sample  quantization.  Expressed  in  decibels,  this  expected  loss 

ratio  is 


Loss  •=  -10  log^Q{(K+2-N)/(K+l)} 

-10  log^o{  (S*Mq^S)/(S*M"^S)  )  OO) 

where  the  first  term  is  the  loss  in  decibles,  due  to  the  finicu  sampling 
in  estimating  M,  and  the  second  tern  is  the  loss  due  to  quantization  of 
finite  step  size  A.  Using  (20),  an  estimate  of  the  last  term  can  be 


made  from  the  identity, 


+  0(a^) 


-1  -2 
K  M  -  OM 


where 


a  “  A /6. 


Hence 


*  -1 
S 


*  -1  A  *  -J* 

S  M  S  -  S  h  i 


and 


^  '  t  K  S 


-10  log,.  II  - 


s  y 


(31) 


is  the  S/N  loss  due  to  quantiidti','!  is  the  optimum  weight 

vector  (21).  The  reliability  cf  14,S8  formulas  will  be  checked 

by  simulation  during  the  next 
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APPENDIX  II 


FORTRAN  PROGRAM  FOR  SIMULATION 
OF  ADAPTIVE  AMTI  WITH  NON-LINEARITIES 


n  ca 
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program  ADAPTM8 
C 

C  MOD  INCLUDES  SAMPLE  MATRIX  INVERSION  SMi, INVERSE  MATRIX  UPDATE 

c  iMu»  and  applebaum  loops 

C  NPaNUMBER  OE  PULSES 

c  nel»number  of  antenna  array  elements 

C  NSCpNUMBER  of  clutter  SCATTERERS  per  range  ring 

C  ALL  distances  ARE  MEASURED  IN  WAVELENGTHS 
C  DpDISTANCE  radar  platform  moves  between  PULSES 
C  ELSPACsELEMEnT  spacing 

C  THETAOalNITlAL  SCAN  ANGLE  IN  DEGREES# ME ASURED  From  GROUND  TRACK 
C  THDOTOaSCAN  RATE  IN  DEGREES  PER  PULSp  REPETITION  PERIOD 
C  RECN*RATIO  OF  RECEIVER  NOISE  TO  CLUTTER  POWERS 
C  SIGSAMPLITUDC  OF  STEERING  SIGNALS  IN  LOOPS 

C  XRaANY  RANDOM  NUMBFK  TO  INITIALIZE  THE  CLUTTER  AND  NOISE  GENERATORS 
C  ULPAT  controls  THE  ELEMENT  PATTERN*** lELPATsi  FOP  ELEMENT  PA 

C  OTHfKWiSE  THE  element  PATTERN  IS  ISOTROPIC 

C  IEI6«1  NORMAL  coordinates  USED  IN  SIMULATION 
C  NPRSSAMPLES  IN  ONE  GROUP 
C  NBP=NUMB£R  OF  GROUPS 
C  ALsWEIGHTING  of  new  data  for  IHU 
C  DLSG  5  DEL/SG 

ITL  SELECTS  CASE  NSCl  VALUES  OF  DLSG  OR  BITS 

NSC2  VALUES  OF  AK  (AK  *  SG  e  SAT) 

ITL  a  0  LIMITING  ONLY  NSCt*!#  DLSGsO#  NSC2  VAIUFS  OF  AK 
ITL  P  I  quant  only  NSCi  VALUES  OF  DLSG#  NSC281,  ak  ,GE,  ? 

ITL  «  %  OUANT  and  limiting  NSCI  VALUES  OF  NUMBER  OF  HITS-AND  NSC2 

values  of  ak  for  EACH  ABIT  VALUE 

call  main 
stop 

END 


SUBROUTINE  MAIN 

COMPLEX  W(T;0)  #  WOI  iO)  #SCAT{60)  #SST(30)  #DCT(30) 

COMPLEX  CAIN(30#60)#V5{6#3?0)#  SC { 30# 3P) # CM ( 3o# 30) # V (8, 320) 
COMPLEX  CN(3o#30)#  AC30#30,)  #B(30»  30)  #CO(30f  30)#CC#SS 
real  LSS#  LTOT 


DIMENSION  EIG(30)#SE1G(30)#CAS1 (20) # CAS2{20) 

DIMENSION  Y1 I (500)» YI 2(500 )#Y2I (300), Y22(500),V31 (SOO) , Y32(S00) 

OIMENSIUNXJ (5OO),SCL(5O0)#X2(iO0)#SCS(lO0)»X3(lOO)#SCi(  100) 

C0MM0N/PARAM/NP#NIL#NPL#NSC#S1G 

COMMON/ARRAY/CAlN# VS#SST 

COMMON/B/WjKO 

COMMON/M/SC 


COMHON/N/CM 


COMMfJN/P/Xi,X2#X3/8CL,SCS;5CI#n,l2,r3#INV#LP#ISP#NRS 

COMMON/O/CO 


COMMUN/R/DCT,ThdCT# AOCT# JDCT 
COMHON/NS/ICX# ICY# ICT 
CUmmOn/Y/YU#  Yi2,  Y2j#y:2»Y31#Y32 

NAMEL IST/INPUT/  NP#NEL#NSC#NRuNG,NSIC#NPR# Inv,lp#isp#nrs# 
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lERR#GLOP,AL,NBP,Df8IO,XR,lELPATiElSPAC,THtTAD,THDDTOrRECN, 
2IEIG, IPUOT/ASIG, aDCTiTHOCTD, JUCT»KDCTf AK»NBIT#DEL#SAT»DLSG 
3,TLDB,NHN5,NCS1,nCS2,1TL,IUT 
REAP(5, INPUT) 

IF{NEI.,EQ,0)  RETURN 

I5AV*IE1G 

Pl«a,*ATAN(l,) 

RADePI/180, 

XRbRANF(XR) 
nRunssnspanpr 
I  fORMAT(lOF5,n 
THDOTpTMDOTD^RAD 
THETAaTHETAO*RAO 
THDCTsTHDCTD*RAO 
NPL«NP*NEt 
WRITF(6»INPUT) 

c  ideal  performance  at  angle  theta 

NRsl 

CALL  P..ASE  fD#CONST»SCZERn,  IELPAT»ELSPAC, THETA,  YHDOT# 

WRITE(6,5^2)  CONST, gCZERD 
SRECNeRECN*NfcL*f  3/2.XSC2ERO*SI6***2 
VRECNsSGRT (SRECN) 

C  IDEAL  COV«  MATRIX 

call  CVM(NpL, sc, NSC, SRECN) 

TAU2GL0P*HEAL(CM(I,l))ANPL/2,/rRR 
PRINT  500, TAU 

PRINT  500, (  SC(i;i)  ,l3l,NPL) 

CCbO. 

DO  12  1=1, NPL 

12  cc*cc+sc(i,n 

CC*CC/NPL/2, 

SGs5QRT(CAHS(CC) ) 

PRINT  530, SC 
PRINT  510 
PRINT  502 

call  WTS(W, SC, const, 0,0) 

C  SC  =  INVERSE  IDEAL  MATRIX 

NR8=«SI3«1 

call  SIGCN(W,5ST, NFL, CONST) 

NRSaS 

call  51GCNCw,S8T;nPL, CONST) 

SUSCKl) 

SCI(1)b10,*ALOG10(SCI(1)) 

SSGsSCKl) 

SCS(l)aSSG 
C  DIAGONAL  OF  SC 

PRINT  510 

PRINT  500, (SC(I,I),IP1,NPL) 

'■  CUMP  vs  FOR  EACH  SAMPLE 
AD*AOCT*NPR»nBP 
NbPR  8  NPK*NBP 
DU  31  J  s  1,NPPR 


0,NR) 
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CALL  SC6EN(NSC,SCAT,XR) 

DO  31  LsliNPL 
CC5O, 

TF(JtEQ.NSlC)  CC  9  A8IG*CONJG(SST(L) ) 
iFf J.EO.JDCT)  CC  *  CC+AD*DCT(L) 

DO  32  KslfNSC 
CC»CC4'SCAT(K)*GAIN{L»K) 

32  CONTINUE 
31  VS(L»J)=CC 

IF(HECN,6T,0,}  CALL  NA00£J.N8PRr VRECN,XR) 
DO  37  J=1»NBPR 
DO  37  L=l»NPL 
37  v£L/JI=VSCLf J) 

CStsO, 

DO  20  Jsl/NBPR 
DO  20  LsUNPL 

80  CST=CST+RLAL(V  (L#J)*C0NJG(V  CLrJ))) 
CST=CST/2, 

8t*SQRTCCST/(N8PR*NPL)) 

PRINT  530, SG 
5^1  *=  0, 

CC  3  0, 

DO  10  I  s  1,NPl 
88  9  SStN(I)*Cf|NJG(w(l)) 
lO  CC  a  CCfSSTCn*CnNJG(W(n) 

ACaRFAL(CC/SS) 

XXsSl/CONST 
PRINT  500, AC, XX 


Ife  CONTINUE 

n6PRkNPH*NBP 

L88=10,*ALOG10((NPR+2,«NPL)/(NPRfJ,) ) 
NBsNBP*  1 
DO  9«  Ial,NB 
9«  SCHDrSSG 

PRINT  509,ITL,ILT 
RtAO(b,l)  (CASl(IJ,lal,NCSl) 

PEAD(5,1)  tCAS2(I),ial,NCS?) 

DO  301  18=1, NCSl 
DLSGsCA81(IS) 

DO  300  IRsi,NC82 
ICXCICY=ICTS0 
PRINT  510 
PRINT  510 

PRINT  SUOfCASl (IS),CAS2(IR) 
AKsCAS2(IR) 

DO  35  M  B  1,NPL 
DO  55  N  B  ) ,nBPK 

35  v8(M,N)  b  vim,n) 

IF  CITL,ED,0)  go  to  22 
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IF(IT1.,EQ,1)  Gf)  TO  41 
ir(ITL,E0,3)G0  TO  3 
TLOBsCASKIS) 

TLOS=lO,**(*CASUIS)/lO,) 

DlSGsSORK  (1,-TL0S)*ACA6,)/SC 
«  CONTINUE 

1F(DLSG,LE.  0,5  GO  TO  <>1 
Nl,s2f*AK  /OLSG+.S 
N8lT3AL0G(NL+2,)/AL0G(a,)  +  ,«> 

AKsDLS6*NL/a, 

GO  TO  5 

5  NBnsCASl(IS)  +  ,5 

NLsa^'^NB  I  T^a 
OLSG=2,*AK/nL 
h  CONTINUE 

DSl3T0LSG*SG)**2>6t 
DO  1«  M5WNPL 
DO  Iti  N=1#NPL 
CC  sCM(M,N) 

IF(H,EQ,N)  CC=CC  tD5l 
H  SC(M,N)aCC 

CA',L  HAriNV(NPL*SC) 

CCsO, 

DO  11  Msl/NPL 
SSsO, 

DO  13  Nsl,NPL 
13  SSaS8fSC(M,N)*SST(N) 

11  CC=CC+CUNJG(SST(M) )*S3 
S5CsR£At(CC)-»CONST 
XXsSSC 

Tl.ni  a  10,*ALOG10(S5C) -SC  un 
LTnT*LSS*TL01 
PRINT  SaSfTLDl.ESSfLTOT 
XX=1 .-OSl/AC 

IPCXX.GT.OJ  D0  a  lQ,tAL0GlQ(A8S(XX) ) 

PRINT  500rXX,DB 
ABlTaNBIT 
ANI  sNt 
23  CONTINUE 
$ATbSG* AK 
PRINT  S26 

PRINT  500»  SAT,AK,TLDB,DLSGf ABIT,ACfCC,S3f ANL 
IK(ILT,£Q,1)  CALI  OUANTI (VS»NL,SAT,NPL»NBPR,ITL) 
IF{ILT,E0,3)  call  QUANT2(V3,NLfSAT,NPL»NBPR,lTL) 
CTcPLOAT(ICT) 

CXsFLOATdCX) 

CYaPLOAKICY) 

PRINT  SH 

PC»(CX+CY)/CT»100i 
PRINT  516, PC 
PR1NT500,CT»CX,CY 
PRINT  510 
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CONTINUE 
IJslSRPlUJsl 
xi(n30$x2(n=os>x3(i)50 
NP|_4s(4*NPI. 

c 

C  BIG  LOOP 

c  simulation  using  random  clutter  samples 

NBOSl 

NCTsO 

JJpO 

DO  87  I  a  IfNBP 
JJ  «  JJ+NPR 


C  8MI 

NRSs2 

I2sl2+1 

X2(12)sJJ 

NUJJ^-NPR+i 

CALL  SAMPLt  (NPL»NPR,I,CONSTrNn 
a  CONTINUE 
87  CONTINUE 
C 

AV  X  0, 

AVS90, 

VA  s  0, 

IZXI2»-1 

DO  28  1  s  2,IZ 
AVSsAVS  +  SCSd) 

B6  AV  =  AV  +  SCKI) 

AV=AV/(N8P«1) 

AV6aAVS/(N8P-l) 

VAsO, iVASxO, 

DO  29  I  B  2,12 
IF(I,EQ,13)  GU  Tn  90 
VAS  =  VAS+(5CS(I)/AVSi-1,)**2 
VA  sVA  ♦(5CI(n/AV  -1,)**2 
90  SCSI  nxlO.AALOGlOCSCSd)) 

89  SCId)xlO,*ALUGlO(SCI(I)) 
VAeVA/(NBP-n 
VAS9VAS/(NBP-n 
AV  «10,*ALOGIO(AV  J-SCKl) 
AVS=10,*ALOGiO(Av8)-SCl(n 
PRINT  50« 

PRINT  500»  (X2d)#l«l»12) 

PRINT  506 

PRINT  5l5fAVS,VAa 

PRINT  500»  (SCSd),I*l,I2) 

AVL30, 

VALPO, 


DO  50  I=3»ia 
50  AVLBAVt  +  SCLd) 

AVL*AVL/(N6P*n 
DO  36  I=3rI2 

VAL«VAL+(5CL(n/AVL*l,)^*a 
36  sCL(n  =  10,*ALOGlO(SCL(I)) 

val»val/(nbp-1) 

AVL  =  10,*ALOGIO(AVL)'SCI(U 

PRINT  515»AVL»VAl 

PRINT  500»  (SCL(n#l«3#l2) 

PRINT  510 
IF(INV,EQ,1) 

IPRINT  507 
iFdNV.EQ.l) 

2PRINT  500»AV,VA 
lFdNV,EQ,l) 

3PRINT  500f  (SCICI)d  =  l#n) 

PRINT  510 

GO  TO  (80f81,6P);  13 

60  Yll (IR)=SSG*AV8 
Yl2(IR)=S5G+AVL 
GO  TO  300 

61  y21 (IR)=SSG+AVS 
Y22(IR)=SSG+AVL 
GO  TO  300 

82  y31 (IR)5SSG+AVS 
Y32(IR)sSSG+AVL 

300  CONTINUE 

301  CONTINUE 

DO  302  I  s  l»NC82 
302  Xl(I)=CAS2(n 

lKIPLOT,EQ,n  CALL  0PUOT(NCS2»SSG) 
yR»RANF(*l.) 

AKsAKS 
read  INI’UT 
1F(NEL,EQ|0)  RETURN 
XRsRANF(XR) 

GO  TO  16 

500  F0RMAT(lX,6E16.:j 

BOl  FORMAT(/»*  NUMBER  OF  SAMPLES  =  *rI5) 

502  FORMAT!/  t*  OPT  performance*) 

503  FORMAT!/;*  XRF  *,F17,10,<^  TAUA  *,E17,10) 

50a  FORMAT!*  X  AXIS*) 

506  FORMAT!*  8AMPLF  MATRIX  INVERS ION ( 8M I , SC: ' * ) 

507  FORMAT!*  INVERSE  MATRIX  UPDATE ! IMV , SC  I ) *) 
606  FORMAT!1X,101A) 

609  FORMAT!/;*  ITLa  *;!«;*  ILTb  A;I«) 

610  FORMAT!/) 

511  FORMAT!/;*  EIGEN  values*) 

512  FORMAT!/;*  IDEAL  COVAR,  MATRIX*) 

$13  FORMAT!/;*  NDRMAl  MATRIX*) 

5!a  FORHAT!lX;*  CTOT  CX  CY*) 


61 


515  FOPMAT(/f 1X»*STMULATI0N  LOSS  «  * , F 1 0 , 3# * VARl ANCL  b  *,F10,3) 

$16  FORMAK/#  lX#*PERCtNT  LIMITING  s  *,F10«3) 

525  F0RMAT(/»1X,*TH,  quant  LOSS  DB  =  *,F10,3»*TH,  SAMPLE  LOSS  s  *,F10, 
n»*T0TAL  TH,  LOSS  »  *,F10»3) 

5?6  FORMATt/#*  SAT  AK  TLDB  PLS6  ABIT  AC  TC  SS«) 

527  FORMAT(1X,10#F12.3) 

630  FORMAT(/f*  SGs  *,2£l6,7) 

S«2  f0RMAT(7H  CUNST»,E16,6# 10X#7HSCZERDs,E16,6//) 

END 


SUBROUTINE  BPL0T(NSP,SSG) 

COMMON  /PARAM/  Np, NtL r NPL # NSC # SIG 

dimension  Yll (500)# Yl2(SO0)#Y21(5OOJ,Y22(bOO),VSl ( 500) , V52 (500) 
OIMENSIONXI  (500),SCL(500)#X2(|00)r5CS(  JOO)rX3(100)f  SCKIOO) 
OlMtNSIDN  AM(500) 

C0MM0N/P/Xl#X2,X3/SCL#SCS,SCl#n»  12#  I3#INV»LP#  ISP#NRS 
r,UMHON/Y/yn#Y12,Y2l#Y22#Y31#Y32 
DO  3  I  =  1»NSP 
3  AM(I)  =  SSG 
NSlFNSP^l 
ns?=nsp^2 
Xl(NSl)s2, 

Xl(NS2)5l  ,0 

Yll  (Nsn  =  io, 

YlPtNSl )siO, 

Y2l(NSl)slO, 

V22(NSl)slO, 

Y31 CNSl)BlO. 

Y32(NS1)=10. 

YU(NS2)s5, 

YIP(NS2)=5, 

Y21(NS2)s5, 

y22(NS2)=5, 

Y3J (NS2)b5. 

Y3?(NS2)s5, 

AM(NS1)=10, 

AM(NS2)e5, 

call  PLDTTSCf7HMALLETT.  U) 

call  AXIS(0,i0,f 7HSAT/SIG#*7»6,#0.»2.#1 ,0#0) 

call  AXIS(0,,0,#12HmTI  GAlN(DU),l2,8,,90,#10,r5,#*l) 

HT*(SSG-Y32(N51))/V32(NS2)+0,1 

CALL  SYMGOL(3,0#HT  , ,l,2JHSTtA0Y  STATE  CAIN(DB)  s#0,»23) 

CALL  NUMeER(5,0#HT  # . 1 # SSG,0, #aHF6. 1 ) 

CALL  LINE(XWAM,n8P»  1#0,0) 

call  line(xi,yji,nsp, 1,1,1) 

call  LINE(Xl,Yt2,NSP, I, 1,2) 

CALL  LINE(Xl,Y?i,NSP, 1,1,1) 

CALL  L1NE(X1,Y?2,NSP, 1, 1,2) 

CALL  LINE(Xl,Y31,N6P#l#l»n 
call  LINE(x1,Y52,N8P, 1,1,2) 

CALL  PLOT(10,,0,,»3) 


RETURN 

END 


SUBROUTINE  PHASE(D,CONSTrSC2ERO»  IELPAT»ELSPAC,  THETA  #THDrjT  fNC,NR) 

c  computes  return  prom  each  scatter  to  each  element 

COMMON/P ARAM/NP,NEL,NPt, NSC #SJG 
COMMON/ARRAY/GAIN,  VSiSST 
C0HM0N/K/UCT,THDCT#ADCT» JOCT 

COMPLEX  GAlNl3r,6O),VS{0»32O)rDCT(3O),58T(3O),CP 

Pls«,*ATAN{l,) 

CPsCMPl,X(0,,?,*Pl) 

CZEROsO, 

THM1N=THETA*PI/2, 

THHAXsTHETA+PI/2, 

CS(THMAX»-^HMIN)/nSC 
KKsINT( (THDCT-TriMiNJ/C+l) 

TH0CT=THMIn+(KK«.5)*C 

NSeN5CflNTC(NR*THDOT)/C)  tl 

PST=THETAfNC*THDOT 

OU  10  K=1»NS 

THsTHhIN*CK«*,S)»C 

CCcCDSfTHj 

DO  10  Ms1,nP 

MNMsMh J  fNC 

8StPSI*TH 

lF(ABS(SS),CT,Pi/2.)  GQ  TO  10 
SS»SIN(SS) 

GTsNEL 

lFfSS,EG,0)  GO  TO  80 

GT  =  SlN(PI»ELSPAC*SS*Ntl)/SIN(PI*ELSPAC<^SS) 

08  CONTINUE 

IF  (lELPAT.NE.n  GO  TO  30 

GAm=,207*COS(PSI*TH)»,256 

ELPsSIN(21,*Pl*GAM)/SIN(PI*GAM) 

GTcGT*ELP**2 
50  CONTINUE 

IF{M,E0,1)  CZEROsCZERO+  GT*Afl 

00  10  Ns1,nEL 

LafMrl)*NEL+N 

GAIN(U  fK)=GT*CEXP(-CP*((N»(NEL+n/2,USS*ElSPAC»MNM  *2,*D*CC)) 
IF(K,E0,KKJ  0CT(L)*GAIN(L,K) 

10  CONTINUE 

SCZER0sNELa*«*SIG**2^CZER0 
CUNST=NEL**2/SCZFR0 
IF(UDCT,EG,0)  GO  T031 
CCrCnS(THl)CT) 

S5sPSI».THUCT 
G  T  s  0  I 

IF fABS(SS),GT, PI/2,)  GO  TO  91 
8S»SIN(SS) 

gt«nel 
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IF(SS,EQ,0)  60  Tn  ^1 

6T*SlN(PI*ELSPAg*5S*NEL)/SIN(PI*ELSPAC*6S) 

91  CONTINUE 
lF(IELPAT,NE,n  GO  TO  51 
THsTHOCT 

GAMS, 207*COS(PS1*TH)., 258 
ELPsSIN(2l,*PI*6AM)/SIN{PI*6AM) 

GT*6T*EI.P**2 
51  CONTINUE 

DO  20  MslrNP 
L0(M»1 )*NEt*l 
MNMsNC*M*2 

SST(L  )cSI6ACEXP(*CP*(  MNM*2,#0*COS (PSl  )  +(M.i)/2,)) 

DO  20  NslfNEL 
LLsL+N 

IFIJOCT.EQ.O)  GO  TO  92 
ISL+N*J 

DCTCI)a  GTACEXP(»CP*(CN*»{N£tin/2,)ASS*EL8PAC*MNM  *2,«»D*CC 

92  CONTINUE 

20  SS7(LL  )=5ST(L) 
return 

END 


SUBROUTINE  NADD(J,NPR,  RECN^XR) 

C  RECi  NOISE  GEN, 

COMPLEX  SST(30),QAIN(30,60), VS(8,320),CP 
common/param/np,nEL,npl,nsc,sig 
COMMON/ARKAY/GAIn# VSiSST 
CPBCMPLX(0,,2,A3,ni5V265) 

00  10  J=l»NPR 
DO  10  MslrNPL 
XR?RANF(0,) 

AMP  =  RECN»SQRT(*.  ALOG(XR)) 

XRsKANF (0,) 

10  VS(M,  J)  =  V5(M,  J)  +  AMP<tCEXP(CPAXR) 

return 

end 


SUBROUTINE  RTS(W, CM, CONST, Ji,J2) 

C  OPTIMUM  WEIGHTS 

COMPLEX  SST(30),w(30),CFPKEC,C,S 
COMPLEX  CM(30,305,GA1N(30,60),VS(6,520) 

COMMON/P ARAM/NP,nEL,NPL,  NSC, SXG 

COMMON/ ARK AY /gain, VS, SST 

S*0, 

HMATaNpANEL 

IF(J1,EQ,C'  call  MATINV(NMAT»CH) 
lFrj2.NE,2)  PRINT  91 

91  FORMAT  (/7X,6MPULSE  ,  9X  ,  7HELEMENT,  6X ,  iiHREAL ,  5X,  9HI MAGINARY,  SX , 
i9HAMPLnuOE,5X,5HPHA3E//) 
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DO  VO  Mslfhjp 
J8(M-1)*n£L 
DO  90  NalfNtL 
MHbN+J 
Coo, 

DO  31  Hl=l,NP 
ie(Ml*l)*NEL 
DO  31  Nlal,NEL 
NN^Nltl 

CsC  +  SST(NN)<(CM(MM/NN) 

3l  CONTINUE 

SsSfC*CONjG(SST(MM  )) 

W(MM  )sC 

lF(ja.EO,2)  GO  TO  90 

AMpsCABSCC) 

PMASATAN2(AIHAG(C)  ,REAL(C)) 

WRlT£(6r50)  M,NrC  »AMp,PhA 

90  continue 

IF(J2,N£,0)  RETURN 
5(1  FORMAT(2ni»5X,aF20,9,2F12,7) 

SRcCABSCS) 

IF(SR,LT,10,E"10)  GO  TO  2 

SSRsSH 

SK*lO,<>ALOGlO(SR) 

CSTal0,*AUOG10(CnNST) 

GO  TO  3 
2  CONTINUE 
SRKsO, 

3  SRXsSR+CST 

PRINT  52»SR,SRX 

52  FORMAT(//5X,7H^s(Db)»,Gl5,B»^^lHS/CACONSTlDB)e,G15,8) 
RETURN 
FND 


subroutine  SCGEN(NSC#8CAT,XR) 

c  NU!SE  gen,  for  each  scatterer 

CO  °LEX  CP,SCAT(60) 
C/^sCMPLX(Ci,,2.*3,Hl592bS) 

UO  10  Nsl»NSC 
AMPsRANF (0, ) 

XHsRANF (0, ) 

10  SCAT(N)«SQRT(-*  ALOG(AHP))aCEXP(CP*XR) 
RETURN 

end 


SUBROUTINE  SIGCN(W/SST»NPL*CONST) 

C  COMPUTE  SIG/  GLUT  USING  WT8  ON  IDEAL  MATRIX 
COMMON/N/CM 

COHMON/P/Xl  ,X2,X3#SCL,SCS,StI,  n,  12,  Hr  INV,LP»  ISP»NRS 
OlMENSIUNXH500),SCL(500)fX2(l00),SCS(l  00  ),X3  (  1  00  ),SCI(l00) 
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COMPLEX  CM(50r30)f l^{30),SSr(30)/CC,SS,CB 
SSsO, 

CCbO, 

DO  71  MsWNPL 
SSsSSiWlM)*CONJCi(5ST(M)) 

C05O, 

DO  70  NslfNPL 

70  CBaCB^  W(N)<(CM(M,N) 

71  CCeCCfCl3*CONJG(W(M)) 

SSRsCABS(SS)**2 

CCRsCABSCCC) 

SCn=SSR/CCR*CONST 

SCRsSCO 

iFCNRS.cQ.l)  SCL(I1)»SCR 
IF(NR3,LQ,3)  SCI(13)=SCR 
IF(NRS,£Q.2)  SC5(I2)sSCK 
IF(NRS,NE,<4)  GO  TO  2 
8COBslO.*ALUG10CsCtn 
SSRalO,*ALOGlOCSSH) 

CCR=lO,*ALOGlO(CCf<) 

WRITE(6#80)  SSR»CCR»«CDB 

80  f0RMAT(//8H  SIGNAL=^FlO,5»aX»lflMCLUTT£R  POWER* , F 1 0 , 5# «X / 8HS/C ( DD ) 
liMO.5) 
a  CONTINUE 

return 

END 


SUBROUTINE  CVH(NpL/SC,NSC,SRECN) 

C  FORMS  STEADY  STATE  COV,  MATRIX 
COMPLEX  5bT(30)#DCT(30)fCC 

COMPLEX  CM(30f 30),SC( JO,30),6AIN(iO,60)»VS(B,320) 

COMMON/ARRAY/GAlN# VS»SST 

COMMON/N/tM 

COMMON/W/DCT,THDCT» AOCT,  JDCT 

aradct*adi;t 

DO  10  MsIiNPL 
00  10  NBMfNPL 
c  C  »  0 

lF(JDCT,Nt,0)  CCrCC+A*DCT(N)*CONJG(DCT(M)) 

00  21  JslfNSC 

81  CCsCC  +  GAlNCN,  J)*CUNJG(r.AlN(M,  J)) 

IF{M,E0,N)  CC=CC*SRECN 

CH(M,N)=CC$CM(N,M)=C0NJG(CC) 

8C(M,N)=CCSSC(N,M)*C0NJ6(CC) 

10  CONTINUE 
RETURN 
END 


SUBROUTINE  sample (NPL,NpR, I, CONST, Nl) 
COMPLEX  w(30),WO(30),5ST(iO),CC 


H'ffi  M— ’ll  ■ 
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COMPLfSf  COf30,S0)rCN(50,50)#CZ(30#S0)#GAIN(30#60)#VS(6,3^0) 
COMMON/Q/CO 
COHMON/B/W,Wn 
COMMON/ARKAY/GAiNf V3#SST 
C  COMPUTE  sample  COV,  MATRIX 


wasNJ+NPR*! 

Ar  1  /I 

1F(iInE,1)  go  to  fi 

ASl  , 

00  B  MbWNPL 
00  8  NsliNPL 
8  CN(M#N)sOt 

a  CONTINUE 

1F(A,GT,1,)  Asi, 

DO  ll  MslfNPL 
OO  11  NsM,NPL 
CCsO, 

DO  JO  K?N1,N2 

CC»CCtCDN.Jr,(VS(M,K))AVS(N,K) 

50  CONTINUE 

CN(M,N)aCC/NPR 

CZ(H#N)sCN(M,N) 

CN(N#M)sCONJC(CN(M,N5 ) 

CZ(NiM)sCN(N*M) 

11  CONTINUE 

501  FORMATdX#*  sample  INVERSION*) 
call  H78(m,CN, CONST, 0,2) 
call  81GC2{w,S3T,NPL,CONST,I,CZ) 

J  CONTINUE 

RETURN 
END 


SUBROUTINE  SIGC2(W,SST,NPL,CnNST, :,UM) 

COMMON/P/X 1, X2,X 5, SCLrSCS, SC  I, II, 12, 13,1 NV,LP, ISP, NRS 
DIMENSIONXI (500),SCL(5C0),X2(100),SCS( 100),XJ(100),SCI{10C) 
COMPLEX  UM(5C,  50  5  ,^^(30)  ,  SST(50),CC,SS,C8,US( 30,30) 

SSsO, 

CCsO, 

DO  71  Mci,NPL 
SScSStW(M;*CONJG{SST(M)) 

CB»0, 

00  70  ns1,npl 

10  CBsCB*  W(N)*UM(M,N) 

71  CCcCOCBHnNjGlwfM) ) 

SSRcCAeS(SS)**2 

CCRsCAhS(CC) 

SCn?SSR/CLP*CONST 
SCS(I2)aSLO 
IMI.EQ.I)  CO  TO  5 
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cc«o, 

DO  J  H«l,NPi, 

CBxO, 

DO  2  NsJiNPL 

a  CBsCB4  W(N)*US(M,N) 

I  CC  =  CCfCONJG(  W(Mn*C6 

CCRcCABS(CC) 
SCL(I2)sSSR/CCR*CONST 

5  CONTINUt 

DO  6  Msl,NPL 
DO  6  NeliNPL 

6  llS(M#N)BUM(M,N) 

RfcTURN 

END 


SUBROUTINE  OLIANTl  (V,NL,5AT,NPL,NPR,  ITL) 

C  quantization  SUBi^OUTInE  for  A/D  CONVERTER  OF  FIGURE  3 
COMPLEX  V(8#320) 

CQMHON/NS/ICXiICY# ICT 
DO  10  M  s  i,NPL 
DU  10  N  s  l,NPR 
lCTsICT+1 
yXsREAL (V(M,N)  ) 
yAsABStXX) 

YYsAIMAG( V(M,N) ) 

VA*ABS(YY) 

IKXA.LT.SAT)  GO  TO  9 

yX*SAT<«XX/XA 

lCXalCX+1 

9  1F(YA,LT,SAT)  GO  TO  5 
YYaSAT*YY/YA 
icyslCY^i 

S  iFdTL.EQ.O)  GO  TO  10 

2X  B  (XX/2,/SAT  ♦  .MaCnL  )^,5 
ZY  «  (YY/2,/SAT  ♦  .5)*(NL  )+,5 

yX  s  AINT(ZX)/(NL  )*2*SAT*SAT 
YY  a  AINT(ZY)/(NI.  )*2*SAT*SAT 
10  V(«f N)BGMPLX(XXr VY) 

return 

END 


SUBROUTINE  OUANT?(V,NLiSAT»NPL»NPR» ITL) 

C  quantization  subroutine  fur  A/D  CONVERTER  OF  FIGURE  3 
COMPLEX  V(6r320) 

CUmMON/NS/ ICX,  ICv»  ICT 
DO  10  M  a  1,NPl 
DO  10  N  s  1,NPR 
AMPsCABSC V(M,N) ) 
iCTf ICTtl 

IF( AMP.LT.SAT)  Go  iU  5 
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V(M#N)aV(M,N)«f8AT/AMP 

iCXsICX+J 

If-'(ITL,EO,0)  GO  TO  10 
XXaREA|,(V(M,N)) 
YT^AlMA&CVtMiN)) 

ZX  s  CxX/2,/SAT  4 
ZY  P  (Yy/2,/SAT  4 
XX  «  AIMC2X)/fNL 
YY  a  AINT(ZY)/(NL 
JO  V(M,N)sCMPlX(XX,YY) 
RETURN 

end 


»5)*(NL  )4,5 

»5)4fNL  )4,S 

)*2*SAT*SAT 
)^2*SAT«SAT 


SUBROUTINE  MATINV(N#A) 
COMPLEX  A(30,30) 

00  11  Nlsl,N 
DO  12  J?lrN 
1E(J,E0,N1)  go  To  12 
A(N1>J)sA(N1»J)/a(N1#N1) 

12  continue 

DO  IS  1=1, N 

ifci,eo,ni)  go  To  is 

DO  16  J  =  t , N 
lEfJ,Ea,Nl)  GO  TO  16 

A(I#J)sA(J,J)*a(i  ni)*A(N1,J) 

U  CONTINUE 

A(I,Nns,iA(l,Nl)/A(Ni,Nl) 

IS  CONTINUE 

A(N1,N1)s1,/a(n1,N1) 

11  CONTINUE 

return 

END 
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